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I . SUMMARY 

Work  performed  on  Contract  No.  F8606-75-C-0045  has 
been  reported  in  detail  in  a series  of  eleven  topical  and 
quarterly  technical  reports.  This  final  report  summarizes 
the  material  covered  in  each  of  the  technical  reports  and 
discusses  the  conclusions  obtained.  The  primary  objective 
of  the  program  is  to  develop  methods  for  estimating  the 
yield  of  underground  nuclear  explosions.  The  topics  addressed 
include  the  modeling  of  both  single  and  multiple  explosions, 
propagation  of  the  resultant  stress  waves  through  realistic 
earth  structures,  and  prediction  of  short-  and  long-period 
explosion  seismograms  recorded  at  teleseismically  located 
receivers.  The  results  of  these  investigations  provide  a 
theoretical  framework  for  expressing  uncertainties  in  explo- 
sion yield  estimates  in  terms  of  uncertainties  in  the  near 
source  material  properties,  local  source  and  receiver 
crustal  structures,  and  the  upper  mantle  structure  of  the 
earth. 
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II.  INTRODUCTION 

This  final  report  summarizes  work  performed  under 
Contract  Number  F8606-75-C-0045,  entitled,  "Improved  Yield 
Determination  and  Event  Identification  Research",  by 
Systems,  Science  and  Software  (S 3 ) in  La  Jolla,  California. 
The  program,  which  was  conducted  during  the  period  from  1 
May  1975  through  30  September  1976,  consisted  of  the  follow- 
ing major  tasks: 

• A combined  theoretical/empirical  evaluation  of 
the  effects  of  variations  in  near  source  material 
properties  and  emplacement  parameters  on  explo- 
sion generated  body  and  surface  waves  used  for 
establishing  magnitude-yield  relationships. 

• Development  of  procedures  for  predicting  tele- 
seismic  signals  from  explosions  and  estimating 
uncertainties  in  seismically  determined  yields 
in  terms  of  uncertainties  in  gross  earth  struc- 
ture, and  local  source  and  receiver  structure. 

Detailed  results  obtained  for  these  tasks  have  been 
presented  in  a series  of  Technical  Reports.  These  results 
are  summarized  in  the  following  section  of  this  final  re- 
port. All  reports  prepared  under  this  contract  are  listed 
in  the  Appendix. 
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III.  SUMMARY  OF  TECHNICAL  RESULTS 

3.1  TECHNICAL  REPORT  NO.  1 

"Seismic  Coupling  from  a Nuclear  Explosions  The  De- 
pendence of  the  Reduced  Displacement  Potential  on  the  Non- 
linear Behavior  of  the  Near  Source  Rock  Environment,"  Cherry, 
et  al_j_,  September,  1975. 

3.1.1  Introduction 

This  report  describes  the  results  of  a theoretical 
study  directed  at  the  determination  of  the  sensitivity  of  the 
equivalent  elastic  source  from  a nuclear  explosion  to  the  non- 
linear behavior  of  the  near  source  rock  environment.  The  re- 
sults of  the  parameter  study  provide  valuable  insight  concern- 
ing the  importance  of  various  material  parameters  for  accurate 
prediction  of  seismic  coupling.  The  results  also  axd  in  the 
establishment  of  guidelines  for  the  collection  of  geophysical 
data  in  the  near  source  region. 

3.1.2  Theoretical  Considerations 

A relationship  between  the  teleseismically  determined 
body  wave  magnitude,  m^,  and  the  elastic  properties  of  the 
near  source  material  is  given  as 

m^  'v  log  [a  Y (®)  ] (3.1) 

where  'J' (»)  is  the  steady  state  value  of  the  reduced  displace- 
ment potential  and  a is  the  near  source  compressional  wave 
speed.  This  relationship  was  derived  by  considering  the 
impedence  contrast  between  the  low  velocity  source  region  and 
the  basement  rock.  Not  included  in  this  relationship  are  the 
effects  of  pP  or  any  of  the  details  of  the  geophysical  proper- 
ties of  the  crustal  layers. 
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If  and  H^f®)  are  two  values  of  ¥ («)  cor- 

responding to  materials  i and  k,  then  the  change  in  tele- 
seismic  magnitude  is  given  by 


x K 
Am  = m — m 


aL  ¥.(«) • 


(3.2) 


Equation  (3.2)  gives  the  scaling  law  for  teleseismic  magni- 
tudes in  terms  of  a ¥ (®) . This  equation  is  always  valid  as 
long  as  the  spectrum  of  the  equivalent  source  is  flat  within 
the  teleseismic  frequency  band  and  path  effects  associated 
with  earth  structure  are  invariant  between  events.  Finally, 
uncertainties  in  yield  estimates  based  on  experimentally 
determined  seismic  magnitudes  may  be  related  to  uncertainties 
in  near  source  material  properties  via  direct  application  of 
Eq.  (3.2). 


3.1.3  Results  of  Parameter  Study 

The  results  of  the  sensitivity  study  are  summarized 
in  Table  3.1.  Some  general  remarks  are  in  order  before  dis- 
cussing the  effect  of  individual  parameter  variations.  Note 

that  the  ambient  density,  p , and  the  bulk  modulus,  k , are 

o o 

included  in  the  table  although  they  are  not  independent  para- 
meters in  the  context  of  this  study;  i.e.,  their  values  are 
determined  by  the  independent  parameters  appearing  in  the 
first  ten  columns.  Calculation  No.  18  is  selected  as  the 
reference  for  body  and  surface  wave  magnitude  variations  in 
order  to  avoid  the  appearance  of  minus  signs  in  the  Am^ 
column.  Calculations  marked  with  an  asterisk  in  the  ¥ («) 
column  were  performed  with  a new  source  description  which 
gave  an  RDP  42  percent  higher  than  the  RDP  obtained  in  Cal- 
culation No.  8 using  the  old  source  description;  the  RDP's 
calculated  with  the  new  source  description  were  therefore 
scaled  down  by  a factor  of  0.704  and  the  scaled  values  are 
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reported  in  Table  3.1.  All  calculations  appearing  in 
Table  3.1  were  performed  with  a device  yield  of  20  tons. 

Since  the  elastic  properties  of  the  near  source 
material  appear  in  the  magnitude  relation,  Eq.  (3.2),  via 
the  P wave  velocity,  a,  it  is  obviously  important  to  deter- 
mine their  values  accurately.  As  can  be  seen  from  Table  3.1, 
a positive  error  in  determination  of  a would  cause  body  wave 
magnitude  to  be  over-predicted. 

The  steady  state  value  of  the  RDP,  4* (®) , which  appears 

in  the  magnitude  relation,  Eq.  (3.2),  is  dependent  on  the 

shock  response  of  the  near  source  material.  Referring  to 

Table  3.1,  the  near  source  material  properties  which  have 

the  most  pronounced  effect  on  the  shock  response  are  the 

water  mass  fraction,  f , and  the  air-filled  void  fraction, 

4>  . Positive  errors  in  either  f or  <J>  could  lead  to  sub- 
o Wo 

stantial  under-prediction  of  seismic  magnitudes.  The  varia- 
tion of  seismic  magnitude  as  a function  of  $ is  presented 

0 

in  Figure  3.1.  Note  that  the  introduction  of  air-filled 
voids  causes  a drastic  reduction  in  seismic  coupling  and  the 
importance  of  establishing  this  parameter  preshot  is  obvious. 

Seismic  magnitude  is  not  very  sensitive  to  Pe,  the 
elastic  pressure,  and  P , the  crush  pressure,  indicating  that 

C 

details  in  the  porosity  model  are  relatively  unimportant. 

Seismic  magnitudes  are,  however,  very  sensitive  to 

the  parameters  which  describe  the  failure  surface.  If  Y , 

o 

Y . or  V are  varied  such  that  the  material  strength  is  en- 
m m 

hanced  the  coupling  efficiency  of  an  explosive  device  is 
impaired.  Thus,  a positive  error  in  the  material  strength 
would  lead  to  an  under-prediction  of  seismic  magnitude. 

■ 


3.2  TECHNICAL  REPORT  NO.  2 


"Teleseismic  Coupling  from  the  Simultaneous  Detonation 
of  an  Array  of  Nuclear  Explosions,"  Cherry,  et  al . , February, 
1976. 

3.2.1  Introduction 

The  problem  addressed  in  this  report  involves  the  cal- 
culation of  teleseismic  ground  motion  produced  by  three  15  kt 
nuclear  sources  detonated  simultaneously.  The  three  sources 
were  equally  spaced  165  meters  apart  and  were  assumed  to  be 
contained,  i.e.,  coupling  effects  due  to  cratering  were  not 
included  in  the  analysis. 

In  an  earlier  series  of  calculations  reported  by 
Cherry,  et  al . , (1975b),  the  region  of  tension  failure  re- 
sulting from  the  shock  interaction  of  the  three  explosions 
engulfed  most  of  the  calculational  grid,  allowing  only  for 
qualitative  estimates  of  the  effects  on  seismic  coupling. 

An  improved  grid  was  used  for  subsequent  calculations 
described  in  Report  No.  2.  The  discussion  that  follows  is 
based  on  results  from  Report  No.  2. 

3.2.2  Near  Source  Ground  Motion 

The  techniques  used  to  predict  the  teleseismic  ground 
motion  from  our  explosive  array  require  that  the  divergence, 

7 • is,  and  curl,  V x 2,  of  the  displacement  field  be  monitored 
over  a spherical  surface  which  is  outside  the  nonlinear  mate- 
rial response  region.  The  time  histories  of  V • 2 and  7 x ? 
are  then  decomposed  into  a multiple  expansion  in  spherical 
harmonics.  This  expansion  provides  an  equivalent  point  source 
representation  which  may  be  used  to  quantitatively  establish 
the  seismic  coupling  of  the  array. 
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Figure  3.2  shows  a sketch  of  the  explosive  array,  the 
nonlinear  region  and  the  spherical  surface  over  which  V • 2 
and  7 * ? were  monitored.  The  computational  grid  extended 
far  enough  beyond  this  surface  so  that  external  reflections 
would  not  obviate  the  free  field  assumption. 

The  tensile  fracture  pattern  around  two  adjacent 
cavities  6.85  msec  after  detonation  is  shown  in  Figure  3.3. 
Note  the  "pre-splitting"  between  the  cavities  and  the  altera- 
tion of  the  fracture  pattern  from  that  expected  from  a 
spherically  symmetric  explosive  source.  This  last  item  is 
illustrated  dramatically  in  Figure  3.4  which  shows  the  dis- 
tribution of  fractures  in  the  entire  grid  at  11.57  msec. 

The  fracture  pattern  above  and  to  the  left  of  the  outer 
source  is  that  expected  from  a single  explosive.  The  altera- 
tion of  this  fracture  pattern  is  caused  by  interaction  of  the 
stress  fields  from  adjacent  explosives. 

3.2.3  Teleseismic  Ground  Motion 

In  order  to  study  the  far  field  seismic  signature  of 
an  explosive  event  like  that  considered  here,  it  is  neces- 
sary to  merge  the  nonlinear  finite  difference  source  cal- 
culations  presented  in  the  previous  subsection  of  this  report 
with  the  elastic  wave  propagation  methods  of  theoretical 
seismology.  This  is  done  via  an  equivalent  elastic  source 
representation  of  the  multiple  explosion  source. 

The  familiar  reduced  displacement  potential,  RDP, 
representation  of  a spherically  symmetric  explosion  is  an 
elementary  equivalent  elastic  source.  For  more  complex 
sources  a representation  in  terms  of  an  expansion  of  the 
outgoing  displacement  field  in  spherical  harmonics  is  appro- 
priate. The  calculation  of  an  equivalent  elastic  source  in 
this  latter  form  for  a complex  asymmetric  explosion  source 
was  first  accomplished  by  Cherry,  et  al^.  [1975a],  A detailed 
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Figure  3.2. 


Sketch  of  the  explosive  array,  the  resulting 
nonlinear  region  and  the  surface  over  which 
the  properties  of  the  outgoing  displacement 
field  were  monitored. 
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presentation  of  the  theory  and  its  implementation  is  given 
by  Bache,  et  al.  [1975a] . 

In  Figure  3.5  the  amplitude  of  the  transformed  reduced 
£ 

velocity  potential,  |H'(w)|,  from  the  1-D  calculations  is  com- 
pared to  the  equivalent  quantity  for  the  two  multiple  explo- 
sions under  study.  The  linear  superposition  is  done  by 
(1)  scaling  the  source  as  if  it  were  a single  45  kt  explo- 
sion and  (2)  considering  the  event  to  be  the  superposition  of 
three  15  kt  explosions.  Due  to  the  cube  root  scaling  laws, 
the  only  difference  in  the  two  is  that  the  second  representa- 
tion is  somewhat  richer  in  high  frequencies. 

From  Figure  3.5,  we  see  that  the  "equivalent  RDP"  from 
the  spherically  symmetric  portion  of  the  multiple  explosion 
calculation  does  not  differ  substantially  from  that  obtained 
by  a linear  superposition  of  one-dimensional  sources.  It  is 
interesting  to  note  that  the  corner  frequency  for  the  no  ten- 
sile failure  calculation  falls  between  that  for  the  3 * 15  kt 
and  45  kt  superposition  sources.  For  the  second  calculation, 
with  tensile  failure,  the  "equivalent  RDP"  is  very  little 
different  than  that  for  the  45  kt  source. 

3.2.4  Conclusions 

Calculations  performed  to  date  have  produced  no  sig- 
nificant enhancement  of  teleseismic  ground  motion  from  a 
multiple  shot  array.  This  was  true  in  spite  of  the  fact  that 
the  near  source  fracture  pattern  from  the  array  differed 
significantly  from  that  produced  by  a spherically  symmetric 
explosion  source.  Fairly  large  perturbations  in  the  seismic 
radiation  pattern  did  appear  but  at  frequencies  outside  the 
teleseismic  band. 

Enhanced  tension  failure  did  occur  in  the  multiple 
shot  calculation.  However,  for  this  to  seriously  affect 
teleseismic  ground  motion  it  would  require  an  explosive 
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spacing,  yield  and  material  properties  combination  in  which 
fracturing  from  the  array  is  much  greater  than  that  produced 
by  a single  charge  of  equivalent  yield. 

An  array  in  which  the  important  parameters  are  opti- 
mized to  significantly  alter  the  teleseismic  ground  motion 
could  probably  be  constructed.  However,  it  is  much  more 
likely  that  anomalous  observations  of  teleseismic  ground 
motion  from  an  explosive  array  are  due  to  overburden  effects 
and  high  (low)  coupling  material  properties.  If  these  fac- 
tors fail  to  explain  the  anomaly,  shock  interaction  calcula- 
tions should  be  performed  but  should  be  based  on  a good  esti- 
mate of  material  strength  for  the  rock  at  the  test  site. 

3.3  TECHNICAL  REPORT  NO.  3 

"Constitutive  Equations  for  Fluid-Saturated  Porous 
Media",  Garg,  S.  K.,  October,  1975. 

3.3.1  Introduction 

This  report  describes  the  development  of  constitutive 
relations  for  fluid-saturated  porous  media,  suitable  for  in- 
clusion in  standard  hydrodynamic  codes  (e.g.,  CRAM  or  SKIPPER). 
The  theoretical  formulation  is  based  on  models  for  fluid- 
saturated  rock  aggregates  previously  developed  by  Garg  and 
Nur  (1973]  and  Garg,  et  al.  [1975]  . 

In  the  development  of  the  constitutive  model  it  is 
assumed  that  (1)  there  is  no  relative  motion  between  the 
fluid  and  the  solid,  (2)  no  significant  heat  exchange  occurs 
between  the  solid  and  fluid  phases,  and  (3)  porosity  ♦ is  a 
function  of  only  the  mixture  pressure  Pc  and  the  fluid  pres- 
sure Pj.  The  last  two  assumptions  may  not  be  valid  under 
certain  field  conditions  and  can  be  relaxed,  although  at  the 
expense  of  some  simplicity. 
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The  first  assumption  (no  fluid  diffusion)  is  equivalent 
to  requiring  [Garg,  et  al . , 1975]  that 


X >> 


C k p 
8 0. 

<P2  Pr 


where 


X 'v-  length  scale  of  interest, 

y^  'v  fluid  viscosity, 

p^  ^ density  of  porous  media, 

C 'v  speed  of  sound  in  porous  media,  ' 

0 

k ,v  permeability  of  porous  media,  and 
0 ^ porosity. 

It  is  straightforward  to  verify  that  the  above  in- 
equality is  satisfied  for  many  field  situations  involving 
shock/seismic  wave  propagation  in  fluid-saturated  porous 
media.  This  report  concludes  with  an  iterative  procedure 
for  incorporating  this  new  constitutive  model  into  standard 
hydrodynamic  codes. 


3.4  TECHNICAL  REPORT  NO.  4 AND  TECHNICAL  REPORT  NO.  5 

"Prediction  and  Matching  of  Teleseismic  Ground  Motion 
(Body  and  Surface  Waves)  From  the  NTS  MAST  Explosion,"  Barker, 
et  al. , February,  1976. 

"Comparison  of  Theoretical  and  Observed  Body  and  Sur- 
face Waves  for  KASSERI,  an  Explosion  at  NTS,"  Bache,  et  al. , 
May,  1976. 

3.4.1  Introduction 

These  reports  present  the  results  of  a theoretical  cal- 
culation of  the  teleseismic  body  and  surface  waves  for  the 
underground  explosion  MAST  and  KASSERI,  and  a detailed 
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comparison  of  the  synthetic  seismograms  recorded  at  the 
stations  of  the  Special  Data  Collection  System  (SDCS). 

This  study  is  a comprehensive  analysis  of  both  explosions 
using  computer  modeling  of  the  close-in  nonlinear  ground 
motion  produced  by  each  explosion,  propagation  of  the  re- 
sulting seismic  waves  through  realistic  earth  structures 
and  computation  of  synthetic  seismograms  at  designated 
teleseismic  stations.  One  of  the  important  questions  ad- 
dressed by  this  exercise  is:  What  pre-shot  measurements 
of  the  near-source  medium  are  required  in  order  to  predict 
the  amplitude  of  the  seismic  signal  to  within  some  specified 
range? 

3.4.2  Explosion  Source  Functions 

Fundamental  to  the  prediction  of  teleseismic  ground 
motion  from  explosions  is  the  calculation  of  the  explosion 
reduced  displacement  potential  (RDP)  which  represents  the 
source  coupling  into  elastic  waves.  The  coupling  was  com- 
puted using  the  one-dimensional  (spherically  symmetric) 
finite  difference  code,  SKIPPER. 

A series  of  source  calculations  were  made  for  KASSERI 
and  MAST.  The  KASSERI  explosion  was  detonated  in  ash  flow 
tuff  at  Area  20,  Pahute  Mesa.  The  working  point  was  well 
below  the  water  table.  Standard  measurements  for  density, 
grain  density,  overburden  density,  water  content,  saturation 
and  P wave  velocity  were  available. 

No  material  strength  data  were  available  for  KASSERI. 
Therefore,  we  were  forced  to  estimate  the  material  strength 
from  other  information.  A parameter  study  of  the  effect  of 
plausible  strength  variations  on  the  seismic  signal  was 
carried  out  as  part  of  the  KASSERI  investigation. 

The  amplitude  of  the  reduced  velocity  potential, 

| Y (u) | , is  plotted  in  Figure  3.6  for  four  (146-149)  KASSERI 


Frequency  (Hz) 


Figure  3. 


. The  source  function  amplitudes  for  the  tuff 
sources  of  Table  3.1  and  the  rhyolite  source 
used  for  MAST.  The  source  functions  are  all 
scaled  to  1000  kt. 
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sources  and  the  rhyolite  source  used  for  MAST.  The  ¥ (u)  is 
essentially  the  far  field  displacement  spectrum.  In  fact, 

z log  a ¥ (1  Hz) , 

«■ 

M s log  y ¥ (0.05  Hz), 
s 

as  has  been  demonstrated  theoretically  in  a number  of  past 
S3  reports  [e.g.,  Bache,  et  al.,  1975a,  1975b]. 

Comparing  calculations  146  and  147,  we  see  that  the 
low  strength  portion  of  the  strength  curve  significantly 
narrows  the  peaked  portion  of  the  spectrum  while  having  a 
relatively  minor  effect  elsewhere.  Comparing  147,  148  and 
149,  as  the  strength  is  decreased  the  spectrum  gets  larger. 
Also,  the  peak  becomes  narrower  and  moves  to  lower  frequen- 
cies. The  MAST  source  exhibits  behavior  between  that  of 
147  and  148. 

For  MAST  and  KASSERI  the  important  parameters  con- 
trolling the  teleseismic  coupling  are  a = 4.2  km/sec, 
y = 169  kbar  for  MAST  and  a = 3.1  km/sec  and  y = 90  kbar 
for  KASSERI.  These  values  together  with  the  source  functions 
in  Figure  3.6  give  a first  estimate  of  the  relative  size  of 
the  body  and  surface  waves  for  these  two  events.  Detailed 
comparison  of  theoretical  and  observed  seismograms  is  made 
in  the  following  section. 


3.4.3  Teleseismic  Body  Wave  Predictions  for  KASSERI  and  MAST 

The  teleseismic  body  wave  comparisons  are  shown  in 
Figure  3.7.  For  each  of  the  five  SDCS  stations  at  which 
theoretical  seismograms  were  made,  we  show  first  the  com- 
parison of  predicted  and  observed  seismograms  for  MAST  and 
then  for  KASSERI.  The  important  result  is  that  the  agree- 
ment at  each  of  the  individual  stations  is  woll  within  a 
factor  of  two  in  amplitude. 
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Figure  3.7.  (Continued) 
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The  source  crustal  structures  over  the  top  two  kilo- 
meters were  taken  from  the  respective  Containment  Evaluation 
Panel  (CEP)  reports.  At  greater  depths  the  structures  were 
taken  to  be  the  same.  For  the  structure  at  the  receiver  an 
average  crustal  model  having  little  effect  on  the  seismo- 
grams was  used  at  all  the  SDCS  stations.  For  the  upper 
mantle  we  chose  a slightly  modified  version  of  the  Helmberger 
and  Wiggins  [1971]  model  HWNE . 

The  comparison  between  theoretical  and  observed  ampli- 
tudes from  the  seismograms  of  Figure  3.7  is  shown  in  Figures 
3.8  and  3.9.  The  data  in  Figure  3.8  are  for  the  "b"  phase; 
those  in  Figure  3.9  are  for  the  maximum  amplitude  in  the 
first  three  cycles,  the  "d"  phase.  Along  with  the  plotted 
b and  d amplitudes,  the  periods  of  the  cycles  at  which  these 
measurements  were  made  also  appear  on  the  figures. 

From  the  comparison  of  Figures  3.8  and  3.9,  we  can 
draw  the  following  conclusions; 


• The  scaling  of  observed  amplitudes  between  the 
two  events  is  quite  consistent  except  for  RKON. 

We  previously  pointed  out  the  marked  dissimilarity 
in  waveforms  at  RKON  for  these  two  events. 


• The  match  between  theoretical  and  observed  ampli- 
tudes is  about  the  same  for  the  two  events.  The 
exceptional  stations  are  RKON  and  FNWV.  It  will 
be  very  difficult  to  explain  the  RKON  anomaly 
with  present  techniques  since  the  epicentral 
distance  variation  between  events  is  only  0.1 
degrees.  Waveform  similarities  at  WH2YK , CPSO 
and  HNME  discourage  attempts  to  attribute  the 
anomaly  to  near  source  effects.  Even  in  our 
theoretical  mo  el  the  interference  pattern  at 
FNWV  is  enough  to  make  the  b amplitude  behave 
erratically  as  we  see  by  comparing  results  for 
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Figure  3.8.  Comparison  of  theoretical  and  observed  b ampli- 
tudes for  MAST  and  KASSERI  at  five  SDCS  stations 
Seismograms  were  computed  for  two  values  of  T/Q 
as  indicated.  The  KASSERI  source  function  is 
denoted  148  and  the  upper  mantle  model  is  HWNE-3 
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Figure  3.9.  Comparison  of  theoretical  and  observed  d ampli 
tudes  for  MAST  and  KASSERI. 
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the  two  events.  However,  it  is  unlikely  that 


current  modeling  techniques,  no  matter  how  ac- 
curately used,  would  permit  a very  accurate 
duplication  of  observed  records  at  this  station. 

• The  scaling  between  the  two  computed  source 
functions  seems  to  be  approximately  correct. 
Perhaps  the  MAST  source  is  slightly  large 
compared  to  KASSERI. 


• Differences  between  theoretical  and  observed 
amplitudes  are  less  than  a factor  of  two  nearly 
everywhere.  A goodly  portion  of  the  discrepancy 
can  likely  be  attributed  to  inaccuracy  in  the 
upper  mantle  model.  More  extensive  studies  em- 
ploying much  more  data  could  improve  this  model 
and  therefore  the  agreement  at  individual  sta- 
tions. 

• Changing  T/Q  from  0.95  to  1.05  increases  the 
period  of  the  phase  measured  by  s 0.03  - 0.07 
seconds.  The  increase  in  amplitude  is  15  - 
25  percent. 

• The  periods  of  the  b phase  generally  agree  with 
the  observations.  If  anything,  they  may  be  a 
little  too  short  on  the  average.  For  the  d 
phase  the  theoretical  periods  tend  to  be  a bit 
too  long. 

• Clearly,  the  further  we  go  into  the  record,  the 
poorer  the  match  between  theory  and  observations. 
However,  the  further  one  goes  into  the  record, 
the  greater  the  effect  of  factors  other  than  the 
coupling  of  explosion  energy  into  elastic  waves. 
Hence  the  emphasis  on  the  b phase. 
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3.4.4  Surface  Wave  Predictions  for  KASSERI  and  MAST 

The  first  step  in  our  procedure  for  computing  explo- 
sion surface  waves  was  to  select  an  appropriate  average 
crustal  model  for  the  travel  path  between  NTS  and  the  five 
SDCS  stations.  Basically,  we  began  with  published  models 
that  seemed  appropriate,  then  perturbed  them  slightly  to 
bring  the  theoretical  group  velocity  dispersion  curves  and 
travel  times  into  close  agreement  with  the  observed.  The 
observed  curves  were  computed  directly  from  the  MAST  obser- 
vations using  the  Ss  data  analysis  program  MARS  [Savino, 
et  al. , 1975] . The  dispersion  curves  for  the  revised  models 
fit  the  observations  to  within  0.04  km/sec  for  periods  be- 
tween 10  and  30  seconds.  The  models  were  CIT109-A  for  WH2YK, 
MCEV-J  for  RKON  and  MCEV-B  for  CPSO  and  FNWV.  The  published 
models  on  which  they  are  based  are  CIT109  [Archambeau,  et^  al . , 
1969]  and  MCEV  [McEvilly,  1964].  Only  the  model  MCEV-J  dif- 
fers to  any  substantial  degree  from  the  starting  model.  For 
the  attenuation  we  used  the  model  of  Tryggvason  [1965] . 

Comparison  of  the  theoretical  and  observed  Rayleigh 
waves  at  two  (HNME  and  RKON)  of  the  five  SDCS  stations  are 
shown  in  Figure  3.10  for  both  MAST  and  KASSERI.  The  radial 
and  vertical  components  are  both  shown.  The  most  important 
result  is,  as  in  the  case  of  the  body  waves,  that  the  scaling 
between  the  synthetic  KASSERI  and  MAST  seismograms  closely 
matches  the  observed  scaling.  This  is  important  because  the 
body  and  surface  waves  sample  different  portions  of  the 
source  spectrum  and  are  sensitive  to  different  properties 
of  the  source  medium.  The  final  point  to  note  is  that  the 
agreement  between  the  synthetic  and  observed  seismograms  in 
Figure  3.10  at  each  of  the  SDCS  stations  is  well  within  a 
factor  of  two  in  amplitude  and  could  be  improved  with  more 
refined  travel  path  models. 
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3.5  TECHNICAL  REPORT  NO.  6 

"Teleseismic  Verification  of  Data  Exchange  Yields," 
Bache,  et  al. , May,  1976. 

3.5.1  Introduction 

This  report  addresses  the  possibility  of  directly  in- 
verting teleseismic  ground  motion  in  order  to  obtain  explo- 
sion yield.  The  procedure  proposed  for  accomplishing  the 
inversion  is  to  employ  a deterministic  model  which  predicts 
ground  motion  from  a nuclear  explosion.  This  procedure 
makes  use  of  the  information  expected  from  data  exchange 
packages  that  provide  information  on  the  near  explosion 
source  environment.  Additional  parameters  must  be  inferred 
from  other  sources  or  determined  by  experiment  on  geophy- 
sical analogues. 

This  report  discusses  three  major  topics:  Equivalent 
elastic  source  calculations  for  a wide  range  of  rock  types; 
the  teleseismic  amplitude  dependence  on  the  source,  including 
the  effects  of  the  free  surface  reflected  phase,  pP,  and  the 
effect  of  upper  mantle  elastic  properties  on  short  period 
P waves. 

3.5.2  Equivalent  Elastic  Source  Calculations 

Several  source  calculations,  corresponding  to  an  ex- 
plosive device  yield  of  150  kt  at  a depth  of  burial  of  ap- 
proximately 560  meters,  were  carried  out.  The  rock  types 
considered  ranged  from  a saturated  NTS  granite  to  dry  shale 
and  sandstone.  Measured  material  properties  for  each  of  the 
rock  types  were  used  in  the  source  calculations. 

Figure  3.11  shows  the  source  functions.  Since  there 
is  a factor  of  ten  difference  in  amplitude  at  1 Hz  between 
the  maximum  and  minimum  source  functions,  it  is  important 
to  subdivide  these  sources  into  classes.  Table  3.2  provides 
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Figure  3.11.  Source  functions  for  all  rock  types  except  207, 
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TABLE  3.2 


A 


Calculation 

Class 

Description 

Rock  Type 

a 

U 

¥ (1.0) 

207 

I 

Weak,  saturated 

Tuff 

2.4 

40 

19 

201 

II 

Strong,  saturated 

Granite 

5.3 

207 

12.0 

203 

II 

Strong,  saturated 

Rhyolite 

4.2 

170 

13.3 

204 

II 

Strong,  saturated 

Tuff 

3.1 

90 

11.2 

206 

III 

Weak,  dry 

Shale 

3.9 

150 

4.5 

202 

IV 

Strong,  dry 

Granite 

5.3 

207 

2.2 

205 

IV 

Strong,  dry 

Sandstone 

4.3 

204 

2.0 
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this  separation.  This  source  classification  procedure  pro- 
vides a framework  for  using  data  exchange  information.  The 
steps  involved  are  the  following: 

1.  Classify  the  near  source  rock  environment. 

2.  Perturb  the  material  properties  in  the  class  to 
fit  the  given  depth  of  burial,  density,  P wave 
velocity  and,  hopefully,  S wave  velocity. 

3.  Calculate  the  range  of  equivalent  sources  which 
are  possible  in  this  class.  This  range  will  in- 
clude reasonable  variations  in  material  strength 
and  air-filled  porosity  consistent  with  past  ex- 
perience from  similar  rock  types. 

3.5.3  Teleseismic  Amplitude  Dependence  on  the  Source 

The  question  addressed  in  this  section  of  the  report 
is  how  the  source  spectra  translate  into  short  period  seismo- 
gram amplitudes  at  teleseismic  ranges;  that  is,  into  m^  esti- 
mates. In  order  to  answer  this  question,  synthetic  seismo- 
grams were  computed  for  the  range  of  source  functions  repre- 
sented in  Figure  3.11.  The  relative  amplitude  and  frequency 
content  of  these  seismograms  was  dependent  only  on  the 
character  of  the  source  and  source  crustal  structure,  all 
other  propagation  effects  held  constant.  These  seismograms 
were  used  to  examine  the  scaling  of  P wave  amplitudes  as  a 
function  of  source  material  both  with  and  without  the  effects 
of  the  free  surface  included  (Figures  3.12  and  3.13). 

The  important  conclusions  drawn  from  this  study  are: 

1.  Given  the  equivalent  source,  then  teleseismic 
ground  motion  is  directly  proportional  to  a. 

The  single  most  important  near  source  material 
property  is  a. 
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greement  of  measured  synthetic  seismogram  amplitudes  with  theoretical 
oaling  relations.  The  free  surface  effect  is  included. 


2.  If  possible,  use  the  b amplitude  to  obtain 
magnitude  estimates.  For  the  yields  and 
DOB  of  interest,  the  b phase  is  least  contami- 
nated by  pP. 

3.5.4  Effects  of  Upper  Mantle  Properties  on  Short  Period 

Seismograms 

The  distance  range  of  particular  interest  here  is  be- 
tween 20°  and  36°.  In  this  analysis  synthetic  seismograms 
were  computed  for  this  distance  range  with  all  factors  held 
constant  except  the  upper  mantle  velocity  structure.  The 
example  chosen  is  the  MAST  explosion  in  rhyolite. 

The  upper  mantle  models  studied  are  perturbations  of 
models  HWNE  [Helmberger  and  Wiggins,  1971]  and  HWA  [Wiggins 
and  Helmberger,  1973] . The  perturbations  were  not  allowed 
outside  the  envelope  of  all  possible  models  consistent  with 
travel  time,  ray  parameter  and  amplitude  data  for  the  western 
United  States,  as  determined  by  Wiggins  and  Helmberger  [1974]. 
The  changes  from  the  starting  models  (which  are  fairly  dif- 
ferent from  each  other)  were  at  most  0.3  km/sec. 

Synthetic  seismograms  appropriate  for  the  epicentral 
distance  from  NTS  to  the  SDCS  stations  RKON  and  HNME  are 
shown  in  Figures  3.14  and  3.15.  These  seismograms  represent 
the  extremes  of  the  variation  observed  at  these  stations 
when  perturbing  the  models  under  the  constraints  mentioned 
above.  The  signal  at  RKON  is  affected  primarily  by  the 
structure  100  km  above  and  below  the  400  km  discontinuity. 

The  HNME  seismograms  are  primarily  sensitive  to  the  region 
below  the  600  km  discontinuity. 

■t  • 

The  b and  d amplitudes  were  measured  on  each  seismo- 
gram and  are  listed  on  Figures  3.14  and  3.15.  At  the  station 


RKON,  within  the  triplication  range,  the  ratio  of  the  maxi- 
mum to  minimum  amplitudes  is  about  a factor  of  2-2.5,  though 
the  b phase  for  HWA-3B  is  exceptionally  small.  Beyond  the 
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Figure  3.14. 


Variation  in  amplitude  and  waveform  for  three 
upper  mantle  models  for  MAST  at  station  RKON. 
For  all  calculations  T/Q  =1.05  and  a nominal 
LRSM  response  was  used. 
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triplications  at  HNME,  the  variation  is  less,  a factor  of 
1.5-1. 6. 

The  small  variation  in  amplitude  at  the  station  be- 
yond the  triplications  is  reassuring.  It  is,  to  be  sure, 
based  on  a limited  number  of  model  perturbations,  but  many 
of  these  perturbations  were  designed  to  change  the  ampli- 
tude at  the  HNME  range.  Any  model  used  must  be  consistent 
with  the  available  travel  time  data,  etc.,  and  these  data 
severely  constrain  the  velocity  profile.  If  we  assume 
that  parallel  layered  models  adequately  represent  the  upper 
mantle  for  our  purposes,  we  can  be  reasonably  certain  that 
errors  in  our  knowledge  of  the  velocity-depth  profile  lead 
to  errors  in  our  computed  amplitudes  that  are  on  the  order 
of  a factor  of  1.5  to  1.6. 
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IV.  SUMMARY  OF  CURRENT  RESEARCH 

In  this  section  of  the  report  we  will  briefly  sum- 
marize the  results  of  work  in  research  areas  that  were 
initiated  near  the  end  of  the  reporting  period  for  this 
contract.  This  work  has  not  been  described  in  separate 
Technical  Reports. 

4.1  EXPLOSION  GROUND  MOTION  PREDICTIONS 

The  successful  prediction  of  the  teleseismic  body 
wave  signatures  of  explosion  events  involves  the  following: 

1.  Computation  of  the  coupling  of  explosion  energy 
into  elastic  waves. 

2.  Constitutive  modeling  of  the  behavior  of  geo- 
logic materials  over  the  range  of  stresses 
encountered . 

3.  Propagation  of  the  so  rce-generated  elastic 
waves  to  teleseismic  distances. 

In  the  following  subsections  we  will  briefly  outline 
some  results  on  the  prediction  of  ground  motion  from  explo- 
sions in  different  NTS  settings. 

4.1.1  Source  Modeling 

In  Figure  4.1  we  show  the  equivalent  elastic  source 
for  explosions  detonated  in  a variety  of  NTS  emplacement 
materials.  The  source  function  is  Y (u>)  , a Fourier  trans- 
formed quantity  that  is  related  to  the  far-field  displace- 
ment  spectrum,  u(u>),  in  a whole  space  by 

* 4» 

u(“)  - RS  ' 

where  R is  distance  and  a is  P wave  velocity.  The  key  para- 
meters distinguishing  the  calculations  are  listed  in  Table 
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EIGHT  SOURCES 

Yields  for  sailing  frequency  axis: 

242,  246  15  kt  (Rainier  Mesa) 

216,  245  150  kt  (Pahute  Mesa) 

244,  247,  248  100  kt  (YUcca  Flat  Tuff) 

221  40  kt  (Paleozoics) 
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Figure  4.1.  Source  functions  for  the  materials  specified  in 
Table  4.1.  The  amplitude  axis  is  scaled  to  0.02 
kt  while  the  frequency  axis  is  cube-root  scaled 
to  the  yield  indicated  for  each  source. 


R-3038 


4.1.  They  are  the  elastic  properties,  density  (p)  and 

velocities  (a, 8),  the  air-filled  porosity  ($  ),  and  the 

_ 0 

failure  envelope  parameters,  Y , Y and  P . Briefly,  we 

• mm 

have 


- 1 (?)’ 


where  P is  the  pressure  including  the  overburden  and  J'  is 
the  third  deviatoric  stress  invariant.  The  meaning  of  the 
parameters  Y#,  Y^  and  P_n  is  illustrated  in  Figure  4.2. 

In  Figure  4.1  we  see  that  the  source  functions  vary 
over  a considerable  range  of  amplitudes  for  far-field  yield. 

The  body  wave  magnitude  is  related  to  the  source  function 

A 

amplitude  (to  first  order)  by  m^  : o ?(f  ),  where  f s 1 Hz . 

The  m^  is  clearly  dependent  on  local  material  properties  as 
well  as  explosion  yield. 

4.1.2  Theoretical  Seismograms 

The  source  functions  of  Figure  4.1  were  used  to  syn- 
thesize short  period  seismogram  recordings  of  NTS  events  in 
the  appropriate  materials.  The  yield  and  depth  of  burial  of 
these  events  were  also  known.  We  had  previously  developed 
good  velocity  and  Q models  for  the  travel  path  between  NTS 
and  a receiver  station  at  teleseismic  distances  which  ob- 
tained good  recordings  of  these  events. 

The  comparison  between  synthetic  and  actual  seismo- 
grams at  a particular  station  is  shown  in  Figure  4.3  for 
seven  NTS  explosions.  The  synthetic  seismograms  include 
the  interference  effects  of  the  free  surface  reflected 
phase,  pP.  The  agreement  between  the  predicted  (light  line) 
and  observed  (heavy  line)  seismograms  is  quite  striking. 

For  all  seven  explosions,  both  the  absolute  amplitude  and 
frequency  content  of  the  waveforms  have  been  matched,  provid- 
ing confirmation  of  the  deterministic  modeling  approach  adopted. 
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Figure  4.2.  Assumed  relationship  between  the  material  strength 
(Y)  and  the  hydrodynamic  component  of  stress  (P) 
for  small  e/em.  When  P is  used  for  the  abcissa 
the  shock  loading  path  has  a slope  of  l-2a. 
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TABLE  4.1 

MATERIAL  PROPERTIES  USED  IN  SOURCE  CALCULATIONS 

Rainier  Mesa 


p 

a 

B 

0 

Y 

0 

Ym 

Pm 

Source 

1.9 

2.5 

1.3 

1.36 

0.05 

1.45 

2.9 

242 

0.05 

0.35 

2.0 

246 

Pahute  Above 

P 

a 

B 

♦. 

Y 

0 

Ym 

P 

m 

Source 

1.91 

3.28 

1.9 

14.5 

0.115 

3.1 

4.7 

245 

Pahute  Below 

P 

a 

B 

♦. 

Y 

0 

Ym 

P 

m 

Source 

2.19 

3.4 

2.1 

0 

0.16 

4.3 

5.6 

216 

Yucca  Above 

P 

a 

B 

0 

0 

Y 

0 

Y 

m 

P 

m 

Source 

1.78 

1.8 

1.0 

14.0 

0.03 

0.77 

1.54 

244 

0.15 

0.6 

2.0 

247 

Yucca  Below 

P 

a 

B 

*o 

Y 

0 

Ym 

Pm 

Source 

1.86 

2.35 

1.3 

3.0 

0.05 

1.45 

2.9 

248 

Paleozoic 

P 

a 

B 

0 

Y 

« 

Ym 

P 

m 

Source 

2.78 

5.45 

3.1 

0 

0.34 

12 

14.1 

221,  222 
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Pahute  Mesa  Below  the  Water  Table 


Figure  4.3.  Comparison  of  synthetic  and  actual  (heavy  line) 
short-period  vertical  recordings  of  NTS  explo- 
sions at  a range  > 3500  km.  The  delay  time 
between  the  P and  pP  phase  for  the  calculations 
is  indicated  on  each  record. 
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4,2  SIMULATION  AND  DECOMPOSITION  OF  MULTIPLE  EXPLOSIONS 
4.2.1  Introduction 

Recently,  we  began  experimenting  with  computer  simula- 
tion and  decomposition  of  multiple  explosion  scenarios.  The 
objectives  of  this  experiment  are  to  develop  procedures, 
based  on  close-in  seismic  measurements,  for  verifying  the 
number  and  yields  of  individual  explosions  comprising  a 
multiple  event.  We  are  also  interested  in  detecting  explo- 
sions that  are  detonated  concurrently  with  a multiple  event 
but  are  located  outside  the  array. 


4.2.2  Experiment  Design  and  Data 


The  data  selected  for  the  first  multiple  event  scenario 
consisted  of  close-in  seismic  measurements  of  the  single, 
contained,  underground  explosion  MAST  detonated  June  19,  1975, 
in  the  Pahute  Mesa  region  at  NTS.  These  data  were  obtained 
in  digital  format  from  Sandia  Laboratories  in  Albuquerque, 

New  Mexico,  and  include  accelerometer  and  velocity  gauge 
seismograms  recorded  in  the  distance  range  from  ground  zero 
to  about  9 kilometers  from  MAST.  The  digital  data  were 
sampled  at  a rate  of  500  points  per  second. 

The  explosion  and  recording  station  configuration 
used  for  the  simulated  multiple  event  are  shown  in  Figure 
4.4.  The  test  scenario  consisted  of  a linear  array  of  three 
equivalent  yield  explosions,  equally  spaced  (355  m)  and 
detonated  simultaneously.  The  explosion  spacing  of  355  m 
was  taken  as  representative  of  row  cratering  shots  in  the 
yield  range  near  150  kt. 

In  the  numerical  simulations  to  be  discussed  herein, 
actual  vertical  component  velocity  recordings  of  the  MAST 
explosion  were  used.  The  simulated  composite  velocity 
seismograms  were  constructed  at  each  of  the  stations  along 
the  in-line  and  45°  profiles  shown  in  Figure  4.4  by  delaying 
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Figure  4.4.  Explosion  and  station  configuration 
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and  summing  the  actual  MAST  seismograms  recorded  at  the 
(approximate)  corresponding  distances.  The  delays  were 
based  on  the  spacings  between  shots,  the  propagation 
velocity  assumed  and  the  azimuth  to  the  recording  station. 

An  example  of  the  construction  of  a seismogram  for  a 
multiple  event  is  given  in  Figure  4.5.  This  seismogram  is 
for  Station  3 along  the  in-line  profile.  On  the  left-hand 
side  of  Figure  4.5  we  show  the  actual  vertical  component 
velocity  seismogram  recorded  at  Station  3 for  MAST.  The  dis- 
tance range  in  this  case  is  0.912  km.  The  multiple  explo- 
sion time  series  is  shown  on  the  right-hand  side  of  Figure 
4.5.  Comparing  the  original  and  composite  signals  we  see 
that,  except  for  the  change  in  peak-to-peak  signal  amplitude, 
the  superposition  of  the  delayed  signals  from  the  three 
explosions  results  in  only  subtle  changes  in  the  shape  of 
the  waveform.  It  is  quite  likely  that  an  analyst  would  not 
interpret  this  seismogram  as  that  from  a multiple  explosion. 

4.2.3  Experimental  Results  for  the  Profile  of  Stations  In- 

Line  with  the  Explosion  Array 

In  Figure  4.5  we  showed  the  original  velocity  record 
for  Station  3 located  0.912  km  from  MAST.  We  also  showed  the 
composite  or  summed  velocity  record  for  Station  3 on  the 
profile  in-line  with  the  three  shots.  The  delays  are  based 
on  the  spacing  between  shots  (355  m)  and  a velocity  of  3.8 
km/sec.  Hence,  delay  times  of  0.093  seconds  are  appropriate 
(0.355  km/3.8  (km/sec)  * 0.093  sec)  for  all  stations  along 
this  profile. 

A series  of  16  narrow  band  filters  with  center  fre- 
quencies (f  ) ranging  from  25  to  100  Hz  was  applied  to  both 
c 

the  original  and  composite  signal  shown  in  Figure  4.5.  The 
time  series  output  from  these  filters  is  shown  in  Figure  4.6. 
For  the  original  MAST  record  at  Station  3,  three  major 
bursts  of  high  frequency  energy  occur  for  this  band  of 
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Filtered  Signal 
- Original  Signal 


Filtered  Signal 
- Composite  Signal 


Figure  4.6.  Narrow  band  filtered  signals  (as  functions  of 
frequency  and  time)  for  the  original  (a)  and 
composite  (b)  signals  for  Station  3 (Figure 
4.5)  along  the  profile  in-line  with  the  shot 
array. 
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frequencies  (Figure  4.6a).  The  times  are  0.55,  1.2  and  1.55 
seconds.  Since  we  have  not  attempted  to  interpret  these 
records,  we  will  not  attempt  to  discuss  the  origin  of  these 
high  frequency  energy  bursts.  However,  the  important  fact 
remains  that  there  is  a high  signal-to-noise  ratio  (S/N)  at 
these  frequencies.  The  significance  of  this  result  is  shown 
in  Figure  4.6b  when  the  same  set  of  filters  is  applied  to 
the  composite  (three  explosions)  event.  Here  the  three  major 
high  frequency  bursts  of  energy  now  show  triple  peaks  where 
before  there  was  only  one.  Separation  begins  at  about  40 
Hz  and  becomes  increasingly  clear  at  higher  frequencies. 

Using  the  Hilbert  transform  described  in  Savino,  et 
al . [1975]  , envelope  functions  were  constructed  for  each  of 
the  sixteen  filter  outputs.  The  envelope  functions  are  shown 
in  Figure  4.7a  and  b for  the  original  and  composite  signals, 
respectively.  The  important  point  to  note  (Figure  4.7b)  is 
the  clear  separation  of  the  three  arrivals,  corresponding  to 
the  three  explosions  in  the  scenario.  This  separation  is 
even  more  dramatic  in  Figure  4.8b,  where  the  sum  of  the 
envelopes  at  the  different  frequencies  (Figure  4.7b)  is 
plotted  for  the  composite  signal.  A similar  sum  is  plotted 
for  comparison  in  Figure  4.8a  for  the  original  signal. 

The  time  separations  of  the  maximum  power  arrivals 
between  1.0  and  1.4  seconds  in  Figure  4.8b  correspond  as 
well  as  can  be  determined  to  the  time  delays  (0.093  seconds) 
associated  with  the  equal  spacings  of  the  three  simulated 
bombs.  In  addition,  note  that  while  the  amplitude  of  the 
composite  signal  in  Figure  4.5  is  more  than  twice  the  ampli- 
tude of  the  original  signal  in  Figure  4.5,  each  peak  in  the 
composite  envelope  sum  is  nearly  equal  in  amplitude  to  the 
peak  amplitude  of  the  original  unsummed  signal. 

Similar  results  are  obtained  for  Stations  2 and  4 in 
this  profile.  Again,  three  major  high  frequency  energy 
levels  are  observed.  Separation  due  to  the  three  explosions 
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Figure  4.7. 


Envelope  functions  (as  functions  of  frequency 
and  time)  for  narrow  band  filtered  original  (a) 
and  composite  (b)  signals  for  Station  3 (Figure 
4.4)  along  the  profile  in-line  with  the  shot 
array. 
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Figure  4.8.  Sum  of  the  envelopes  shown  in  Figure  4.7  for 

the  original  (a)  and  composite  (b)  signals  for 
Station  3 along  the  profile  in-line  with  the 
shot  array. 
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occurs  at  about  40  Hz  and  becomes  more  distinct  with  increas- 
ing frequency.  Average  time  separation  of  the  triple  peaks 
is  0.093  and  0.093  seconds,  respectively.  The  observed  de- 
lay times  correspond  almost  exactly  to  the  input  delay  times 
(Table  4.2) . 

For  Stations  5,  6 and  7,  high  frequency  noise  dominates 
the  records  and  the  high  frequency  signal  energy,  if  present, 
is  too  low  level  to  obtain  any  definitive  results. 

4.2.4  Experimental  Results  for  the  Profile  of  Stations  45° 

from  the  Explosion  Array 

Composite  or  summed  seismograms  were  synthesized  for 
the  stations  oriented  at  a 45°  azimuth  from  the  explosion 
array.  Delay  times  were  determined  relative  to  the  center 
explosion  position  (number  2 in  Figure  4.4)  and  are  listed 
in  Table  4.2.  The  delay  times  vary  and  are  less  than  those 
for  the  in-line  array.  At  those  Stations  (2-4)  where  noise 
was  not  a problem,  the  input  and  observed  delays  agree  to 
1 msec  or  better. 

4.2.5  Summary  and  Discussion 

Close-in  seismic  data  obtained  from  Sandia  for  the 
underground  nuclear  explosion  MAST  have  been  used  to  simulate 
a simple  multiple  explosion.  This  multiple  explosion  con- 
sisted of  three  equal  sized  (150  kt)  aligned  explosions 
(separated  by  355  m)  detonated  at  the  same  instant.  Two 
profiles  were  simulated,  one  in-line  with  the  explosion  array 
and  the  other  oriented  at  45*  to  the  explosion  array.  Seismic 
records  for  the  multiple  event  were  synthesized  by  simple 
summation  of  one  MAST  seismic  record  to  the  same  record  with 
appropriate  delay  times.  The  determination  of  the  delay 
times  was  based  on  the  orientation  of  the  two  hypothetical 
station  arrays  relative  to  the  explosion  array,  actual  sta- 
tion distances  from  MAST  and  a compressional  velocity  of 
3.8  km/sec. 
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Application  of  a series  of  cusp-shaped  narrow  band 

filters,  with  center  frequencies  (f  ) ranging  from  25  Hz  to 

c 

150  Hz,  to  the  synthesized  multiple  explosion  records 
yielded  the  following  results: 

• Accurate  relative  amplitude  and  time  separation 
between  explosions  was  achieved  at  the  very 
close  stations  and,  in  particular,  at  distances 

of  less  than  1.5  kilometers.  Beyond  1.5  kilometers 
the  signal-to-noise  ratio  is  too  low  in  the  high 
frequency  band  of  analysis  to  obtain  definitive 
results. 

• Separation  between  explosions  was  observed  to 
begin  at  frequencies  3.5  times  the  explosion 
frequency  (inverse  of  the  explosion  delay  times). 

In  general,  separation  of  events  becomes  more 
distinct  with  increasing  frequency. 
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